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Abstract 


29 

30 

31 

32 This paper presents the process used by the Marshall Space Flight Center Natural Environments 

33 Branch (EV44) to quality control (QC) data from the Kennedy Space Center’s 50-MHz Doppler 

34 Radar Wind Profiler for use in vehicle wind loads and steering commands. The database has 

35 been built to mitigate limitations of using the currently archived databases from weather 

36 balloons. The DRWP database contains wind measurements from approximately 2.7-18.6 km 

37 altitude at roughly five minute intervals for the August 1997 to December 2009 period of record, 

38 and the extensive QC process was designed to remove spurious data from various forms of 

39 atmospheric and non-atmospheric artifacts. The QC process is largely based on DRWP 

40 literature, but two new algorithms have been developed to remove data contaminated by 

41 convection and excessive first guess propagations from the Median Filter First Guess Algorithm. 

42 In addition to describing the automated and manual QC process in detail, this paper describes the 

43 extent of the data retained. Roughly 58% of all possible wind observations exist in the database, 

44 with approximately 100 times as many complete profile sets existing relative to the EV44 

45 balloon databases. This increased sample of near-continuous wind profile measurements may 

46 help increase launch availability by reducing the uncertainty of wind changes during launch 

47 countdown. 

48 
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1. Introduction 


50 The National Aeronautics and Space Administration (NASA) has been designing, testing, and 

51 flying manned and unmanned space vehicles since the second half of the twentieth century. 

52 During this time, the Natural Environments Branch (EV44) at the Marshall Space Flight Center 

53 (MSFC) has helped ensure that supported launch vehicles can withstand the effects of the ascent 

54 wind environment over the central Florida region. EV44 accomplishes this task by developing 

55 meteorological databases and providing environmental definitions as inputs to a wide variety of 

56 discipline-specific engineering analyses supporting aerospace vehicle design, test, and operations 

57 activities. 

58 EV44 has traditionally used observations made with weather balloons to statistically represent 

59 the ascent wind environment. However, balloon-based data archives have a number of inherent 

60 limitations for this application. Relative high cost makes high-frequency balloon sampling 

61 impractical, thereby limiting sample sizes and increasing statistical uncertainty. Downstream 

62 balloon drift can lead to potential misrepresentation of the ascent environment in a horizontally 

63 inhomogeneous wind field. Also, balloons have lengthy rise times, which prevent high temporal 

64 resolution assessments of ascent wind environments. This limitation reduces the ability to verify 

65 steering commands late in a launch countdown, possibly leading to larger-than-necessary 

66 persistence margins during rapid wind regime change events such as frontal passages. The latter 

67 problem was recognized by NASA, and in 1988 a 50-MHz Doppler Radar Wind Profiler 

68 (DRWP) was installed at Kennedy Space Center (KSC). During the 1990s, MSFC developed an 

69 algorithm to operationally quality control (QC) wind profiles from the DRWP on day of launch 

70 (DOL) (Schumann et al. 1995), and various analyses using subsets of the DRWP archive were 

71 performed (Wilfong et al. 1993, Merceret 1997, Schumann et al. 1999). However, the Space 
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Shuttle Program (SSP) did not certify the DRWP in the later stages of the program. Such a 

73 certification would have been very expensive and time-consuming since a DRWP database did 

74 not exist during the Shuttle’s design phase. EV44 highly desires to certify the DRWP for use in 

75 the design of future launch vehicles, and one of the first steps in certifying the instrument is to 

76 develop a climatological archive that can be used by the ascent loads and trajectory community. 

77 EV44 has recently developed an extensively QC’d database of DRWP wind profiles covering 

78 the period of record (POR) from August 1997 to December 2009. This new dataset is expected 

79 to be widely applied to engineering analyses for future launch vehicle programs. The database 

80 contains a much larger sample size than those from balloons, and has greater temporal coverage, 

81 providing flexibility in the pre-launch assessments. Given these two qualities, the DRWP 

82 provides the ability to increase confidence in our knowledge of a given wind environment, 

83 reducing unnecessary conservation in both design and operational margins, thereby potentially 

84 increasing launch availability. To develop the database, an extensive QC process was 

85 implemented on wind and spectral output from the DRWP. This paper documents the QC 

86 process and the database’s application. Following a brief description of the DRWP hardware 

87 and data, the QC process is outlined in detail. Then, the resulting database and its current 

88 applications are described followed by a summary and forward work plan. A list of acronyms is 

89 provided in Table 1. 

90 2. DRWP Database Description 

91 Although detailed descriptions of the DRWP’s hardware and data processing algorithm have 

92 been previously documented (Schumann et al. 1999), an overview is provided here for context. 

93 The attributes of the data stream used in the QC process are also described in this section. The 

94 DRWP is located just east of the Shuttle Landing Facility at KSC. Physically, the DRWP 
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95 consists of an irregular octagon-shaped antenna field spread across 15,600 m . Coaxial-collinear 

96 elements are set 1.5 m above the ground plane made of copper wire. These elements are 

97 arranged to send electronic pulses at 49.25 MHz through three beams. One beam is pointed 

98 vertically, and two oblique beams are pointed 15° off zenith at azimuths of 45° and 135° East 

99 from North. All beams have a 3° beam width (Bill Gober, personal communication 2010). An 

100 adjacent equipment trailer houses the radar instrumentation and data acquisition electronics. 

101 Figure 1 provides a photograph and schematic of the DRWP. 

102 To measure wind velocities, the DRWP transmits radio pulses in three beam directions 

103 sequentially and measures the return signal backscattered by temperature and humidity 

104 fluctuations in the atmosphere with scale lengths of about 3 m (Rinehart 2004). The signals are 

105 converted to Doppler power spectra by applying a Fast Fourier Transform over 256 frequency 

106 bins at each range gate. In the beginning of the POR, 1 12 gates were set from 2,01 1-18,661 m 

107 every 150 m. After an upgrade during July- August 2004 (Pinter et al. 2006), 1 1 1 gates were set 

108 from 2,666-18,616 m every 145 m. Also before this upgrade, profiles were generated every five 

109 minutes. After the upgrade, profiles have been generated every three minutes. Once the Doppler 

1 10 spectra are obtained for each beam, horizontal velocities are computed using the Median Filter 

111 First Guess (MFFG) algorithm (Schumann et al. 1999), which is applied to the oblique beams 

1 12 only as the vertical beam is not used to calculate horizontal winds (Wilfong et al. 1993, 

1 13 Schumann et al. 1999). 

1 14 Although the MFFG algorithm has many advantages over the traditional consensus averaging 

115 technique used on other wind profilers (Schumann et al. 1999), the algorithm is not immune to 

1 16 acquiring erroneous data. In many instances, the FG simply associates itself with spectral peaks 

1 17 which do not represent the real wind. If a strong, non- atmospheric signal persists within the FG 
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118 window, then radial velocities would contain contributions from non- atmospheric effects. Other 

1 19 instances of suspect data occur when the signal is too weak to calculate a radial velocity. If the 

120 signal to noise ratio (SNR) is less than -15 dB, then the measured radial velocity is replaced by 

121 the FG. This process of FG “propagating” continues for an individual beam until the SNR 

122 reaches -15 dB. Although the radial velocity profiles are smoothed after at least four first guess 

123 propagations (FGPs), using previously recorded measurements for an extended period can 

124 introduce errors in the radial velocity estimate, especially if a non-atmospheric signal exists near 

125 the signal associated with the wind (Wilfong et al. 1993). Heavy rain could cause spurious wind 

126 output as the DRWP can track the velocity of the raindrops instead of the air. Also, large vertical 

127 motions can violate the assumption of a homogeneous atmosphere used in the horizontal wind 

128 computation. 

129 EV44 has archived DRWP data output which was used to develop the QC’d database. In 

130 addition to the computed horizontal wind speed (m s' 1 ), wind direction (degrees), and altitude 

131 (m), spectral width (SW, m s' 1 ), signal power (dB), noise power (dB), vertical velocity (w, m s' 1 ), 

132 number of FGPs (dimensionless), and an internal shear value (s' 1 ) at each gate and beam have 

133 been archived. Vertical velocity is the radial velocity of the vertical beam, with positive values 

134 indicating downward motion. The number of FGPs is the number of times the FG was 

135 propagated for the particular gate and beam, and the internal shear is the change of the radial 

136 velocities per unit altitude. These fields exist for each day from August 1997 through December 

137 2009, but have not been regularly QC’d. The QC process used to develop this database differs 

138 from the DOL QC process. The DOL process (Schumann et al. 1995) requires near-real time 

139 examination of the Doppler spectra which are not available in the data used for the current QC 

140 process. The current QC process will be described in detail in section 4. 
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141 3. Data Display System 

142 A graphical user interface (GUI) was developed to implement the current QC process. The 

143 GUI contains various functions to perform all QC procedures and save the desired output 

144 (Figure 2). Each DRWP data file contains data for a single year, month, and day. Files are read 

145 sequentially and time-height sections of meteorological parameters are initially examined for 

146 potentially spurious features. In Figure 2, for example, anomalies at very low altitudes appear to 

147 correspond to ground clutter, and streaks of enhanced meridional wind (v) from 10-12 km at 

148 1500-1700 Coordinated Universal Time (UTC) seem incompatible with the surrounding 

149 environment. The automated QC is then run and a new time -height cross section of the given 

150 variable is displayed. To perform the manual QC, a box surrounding the data in question is 

151 drawn and data which are flagged by the threshold are removed. An “undo” function exists to 

152 protect against operator error during the manual QC process. Once the QC process is complete, 

153 the QC’d file and manual QC logs are saved. In addition, comparisons between profiles from 

154 low-resolution (LR) weather balloons and DRWP profiles at a desired time can be performed and 

155 images can be saved as desired. The LR balloon database consists of rawinsondes prior to 

156 October 2002 and the Low Resolution Flight Element (Leahy and Overbey 2004) after October 

157 2002. 

158 4. QC Process 

159 In addition to methodologies documented elsewhere, the process used here to QC the DRWP 

160 contains some unique attributes based on data examination. A number of distinct steps are 

161 performed sequentially with flagged data being removed before the next step is implemented. 

162 Indicators are assigned to each gate to denote if data passed all checks or failed a particular 

163 check. This section describes in detail the QC sequence. 


7 



164 


a. Automated QC Process 


165 The automated QC process contains initial procedures to fill data gaps and screen the vertical 

166 beam. The first step in the automated QC process fills data gaps. If greater than six minutes 

167 existed between adjacent timestamps in the original data, then timestamps were inserted at five 

168 minute intervals in the data gap with all variables in the profile containing the missing data flag. 

169 This procedure ensured a data record at least once every six minutes throughout the POR. The 

170 second step in the automated QC process evaluates vertical beam measurements. Since the 

171 vertical beam is not used to calculate horizontal winds, a valid wind calculation could coincide 

172 with an erroneous vertical beam measurement and be falsely flagged. Data from the vertical 

173 beam were thus removed if a signal or noise power report were missing, if the vertical beam’s 

174 SW exceeded 3.0 m s' 1 , or if a systematic error occurred when the Doppler shift from the vertical 

175 beam was near zero (Merceret and Gober 2009, personal communication). This error appeared 

176 when abnormally high \w\ coincided with relatively low SNR. 

Ill The automated process then performs threshold checks and flags data possibly influenced by 

178 convection. Table 2 presents each check and its threshold in order. The process consists of 

179 checking for unrealistic wind reports and isolated data, perfonning a small median test, and 

180 applying thresholds to oblique beam SW, DRWP internally computed shear, w, FGP, oblique 

181 beam signal power, and convection. The rationale for each check is described below, with 

182 thresholds for checks other than the FGP, small median, convection, and isolated datum 

183 presented in the table. 

184 Several automated checks were based on thresholds, which were derived from Carr et al. 

185 (1995) and Merceret (1997) and modified if necessary based on data examination. After 

186 detecting physically unrealistic wind reports, a check was implemented on the oblique beam SW 
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187 so that the homogeneity assumption used to calculate the winds would not be violated due to 

188 excessive turbulence. The DRWP internal shear is useful for detecting large objects in the air 

189 such as airplanes (Merceret 1997). Over Florida, |w| is generally very small, so any large 

190 perturbation in w indicates some anomaly in the air flow or that the DRWP is measuring the 

191 velocity of raindrops instead of the air. The threshold selected here is more restrictive than that 

192 in Merceret (1997) to flag additional convective situations, especially after August 2004. The 

193 meteorological shear check serves the same purpose as the DRWP shear check but it applies to 

194 the zonal ( u ) and meridional (v) wind components. Missing signal power indicated that the 

195 DRWP did not receive a signal at that gate. Note that no check exists for the oblique beam noise 

196 power. An analysis was perfonned showing that missing noise values corresponded to some 

197 erroneous vertical beam SW and velocity reports. However, no such effect existed when 

198 examining oblique beam SW and radial velocities. 

199 The small median check (Carr et al. 1995) flags observations which significantly differ from 

200 their nearest neighbors. The check compares a wind speed observation at a given time and 

201 altitude to the eight observations surrounding it and was only performed if the wind speed of 

202 interest and at least three neighboring observations existed. Thresholds were applied following 

203 Merceret (1997). Once all automated checks were performed on data for the day, gates with no 

204 surrounding measurements were removed to enhance the continuity of the database. 

205 Two additional QC algorithms were developed specifically for the DRWP database, and 

206 differed significantly from the literature. These checks involved testing for convection and 

207 developing a criterion for the FGP threshold, and are presented in greater detail in the 

208 subsections below. 

209 


9 



210 


1) Convection 


211 

212 The convection algorithm is derived from previous work using the 915-MHz DRWP network 

213 at KSC. Lambert et al. (2003) developed a discriminant function based on w and SNR which 

214 had two classes: convection and no convection, and was effective on the 915-MHz DRWP data. 

215 However, because the DRWP is much less sensitive to rain than the 915-MHz DRWP and w 

216 differs from the boundary layer to the free atmosphere, this discriminant function was 

217 determined not to suit the DRWP QC process. In addition, a given parameter at an individual 

218 gate might have the same output in different situations throughout the year. Therefore, the 

219 convection algorithm’s parameters were derived for each month. 

220 SW and w were used to detennine if convection existed at a particular gate. Figure 3 shows w 

221 and SW for 20 August 2001, a day with typical summertime convective activity over the KSC 

222 region. Rain gauge data were obtained from the TRMM website 

223 (http://tnnm.ksc.nasa.gov/trmm/rain/2001/08/DAILY_RPT_AUG20.HTM). The rain gauge at 

224 the field mill closest to the DRWP recorded 0.28 inches of rain during 0300-0500 UTC and 

225 1.42 inches of rain during 2100-2200 UTC. Concurrent SW and w show clear signs the 

226 convection could be occurring: SW increased from just less than 1.0 m s' 1 to near 2.0 m s' 1 , and 

227 \w\ increases from near 0.0 m s' 1 to approximately 1. 5-2.0 m s' 1 . Other meteorological variables 

228 did not vary as significantly, so they were not used to discriminate between convective and non- 

229 convective cases. 

230 SW and w were used in a supervised classification technique to detennine if convection 

231 existed. First, the classes “convective”, “possibly convective”, and “not convective” were 

232 


chosen to classify the convective environment of each gate. Next, training samples representing 
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instances of each class were selected for each month across the POR. Instances where large SW 
corresponded with large \w\ across extensive vertical regions were selected as “convective”, and 
instances of small SW and \w\ were selected as “not convective”. Training samples for the 
“possibly convective” class were selected to help mitigate a false positive “convective” 
classificaiton. A range of 1,638 to 6,428; 1 1,546 to 32,265; and 7,528 to 19,935 training samples 
existed per month for the “convective”, “possibly convective”, and “not covective” case, 
respectively. These sample sizes are of one to two orders of magnitude over what is considered 
highly desirable in the literature (Richards 1993). 

The training samples were used to develop a descriminant function, which classifies a pixel as 
convective or not convective. The descriminant function is a quadratic surface described by 


\W SW„ SW,. SW, 


]* L ; 4 


w sw„ sw„ sw„ 
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02,2 

02,3 
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04,1 
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where K, L, and Q are coefficients corresponding to the covariance of the training samples for 
each class combination and month. SW e , SW„, and SW„ are SW from the east, north, and 
vertical beams, respectively. The training samples were provided as input to the MATLAB 
discriminant function routine (http:/www. mathworks.com/help/toolbox/stats/classfiy.htm) along 
with the data to be classified (i.e., the DRWP data for the day of interest). The routine returned 
the class of each gate, the coefficients of the discriminant function for each class combination, 
and the posterior probabilities of the gate belonging to its detennined class. If DF were positive 
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253 for the “convective / possibly convective” combination and the “convective / not convective” 

254 combination with a posterior probability of at least 0.95, then the gate was classified as 

255 “convective”. 

256 Plots showing gates flagged by the convection algorithm were examined and data were 

257 removed manually based on the extent of the flagged gates and the characteristics of the 

258 corresponding wind field. In Figure 4, flagged data over extensive vertical regions which 

259 corresponded to anomalies in v around 0400 UTC and from 2030-2200 UTC were considered to 

260 be convective and were thus removed. However, some flagged gates neither span extensive 

261 horizontal or temporal regions, nor seem to correspond to anomalies in the wind field. These 

262 gates were not removed as convection was determined not to affect the continuity of the winds 

263 on a large enough scale. 

264 

265 2) FIRST GUESS PROPAGATION 

266 

267 A unique FGP check was developed to suit the database’s multiple applications and to 

268 determine how propagating the FG velocity affects the wind estimate. Using radial velocity 

269 estimates with one FGP is basically equivalent to using estimates which are five minutes old. In 

270 situations when the wind changes little, an FGP should not significantly affect the wind estimate. 

271 However, in dynamic conditions, even a small FGP could lead to an inaccurate wind estimate. 

272 Previous research has thus used limits on the number of FGPs for their respective analyses. 

273 Merceret (1997) used a limit of six FGPs; however, no rationale was given for this threshold 

274 other than using it to prevent wind estimates greater than 30 minutes old from being incorporated 

275 into the analysis. Schumann et al. (1999) used a limit of two FGPs to relate to the capability to 
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276 distribute DRWP data every 15 minutes to the end user. Using the latter criterion, data with 

277 greater than two FGPs would contain wind estimates at least 15 minutes old, and the end user 

278 would not be provided with a new wind estimate. This rationale implies that the FGP threshold 

279 should be selected based on the end user’s application. However, it was desired to have a single 

280 FGP threshold for the current DRWP database for three reasons: First, the FGP threshold 

281 significantly affects the number of available profiles. Second, the current DRWP database can 

282 be used for applications with varying time separations which are currently unknown. Last, the 

283 exact relationship between number of FG propagations and resulting measurement errors is 

284 unknown. Therefore, the following analysis was perfonned to better understand the effect of FG 

285 propagation on the output of the MFFG algorithm. 

286 Spectra from three days obtained from the DRWP operations and maintenance contractor 

287 during 2009 were examined. Each day represented a weather regime common to eastern Florida: 

288 • 21 August: Light winds with an afternoon thunderstorm 

289 • 21 October: No rain, dynamic day with moderate winds 

290 • 4 December: Strong southwest winds 

291 Control wind components from each day were calculated using FG radial velocities derived from 

292 the Doppler spectra. First, a three-point median filter was applied to spectra with at least two 

293 valid timestamps. Radial velocities were then calculated for each oblique beam, with the 

294 calculated radial velocity profile being the FG for the next radial velocity profile. The first 

295 profile in the database and the first profiles after data gaps were computed. Radial velocities 

296 were replaced by the mean of the radial velocities from the adjacent gates if the shear criteria 

297 were violated (Taylor et al. 1993). Wind components were then computed from the radial 
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298 velocities, and were compared to the wind components in the output data files for accuracy. If 

299 the magnitude of the difference between the calculated wind component and the wind component 

300 in the output file exceeds 2.0 m s' 1 , then the wind value from the output file replaces the 

301 calculated wind from the spectra. 

302 The control wind components were then differenced from wind components calculated after 

303 propagating the FG velocity from each beam. To simulate propagation of the FG velocity, radial 

304 velocities were calculated using the spectra for the current timestamp with the FG radial velocity 

305 from the previous n timestamps, where n was incremented from 1 to 20. Modified wind 

306 components were then calculated using the previous n FG velocities from each beam. FGPs 

307 from each beam were simulated by cycling n for the north beam before incrementing the east 

308 beam (0 East FGPs / 0 North FGPs, 0 East FGPs / 1 North FGP... 20 East FGPs / 19 North 

309 FGPs, 20 East FGPs / 20 North FGPs). A median of 29,554 gates were used, varying from 

310 28,310 to 32,820 gates depending on the FGP combination. The vector differences between the 

311 modified and control winds were then computed at altitudes over 10 km as the FG is propagated 

312 more frequently at higher altitudes. The RMS of the vector changes for each FGP combination 

313 was then plotted versus the number of FGPs from the oblique beams (Figure 5). Warm (cool) 

314 colors represent large (small) RMS differences from the control wind. As expected, differences 

315 increase as the number of FGPs increase from either beam. An RMS vector error of V2.0 m s' 1 

316 was selected as the threshold for this analysis to correspond to the RMS measurement error 

317 specification of 1.0 m s' 1 for each wind component (Pinter et al., 2006). A quadratic fit was then 

318 applied to the maximum number of north beam FGPs which yielded an RMS difference below 

319 the threshold for each east beam FGP. The fit is expressed as 
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320 T = -0.0 10 (FGP e ) 2 - 0.784 (FGP e ) + 20.309 - FGP n 

321 where T is the threshold parameter and FGP„ and FGP e are the FGPs from the north and east 

322 beams, respectively. If T were less than zero then data at the gate were removed. 

323 

324 b. Manual QC Process 

325 

326 Once the automated process was complete, data for each day were manually examined for 

327 temporal and spatial inconsistencies. This process involved examining multiple variables to see 

328 if a spurious output from one variable coincided with that of another variable. If spatial 

329 discontinuities in multiple variables occurred within the same time-height region, then greater 

330 evidence would be presented to remove the data in question. In addition, temporal changes in 

331 wind components at each altitude were examined to detect the edges of radar sidelobes and 

332 ground clutter. Each manual QC was logged for reference. Data judged to be contaminated by 

333 convection or ground clutter were assigned their own QC flags to be tracked separately. On 

334 occasion, extensive time-height regions of data were not flagged by the automated QC process 

335 but needed to be removed. In these cases, entire time-height boxes or profiles were removed 

336 manually. 

337 An example of the manual QC process is presented here. Figure 6 shows before-and-after 

338 images of v on 19 October 2008. The left panel shows v from the original database. Note the 

339 bar-like features which do not compare well with the surrounding environment at approximately 

340 5.0 km during 0100-0300 UTC and at approximately 4.0 km during 0300-0400 UTC. These 

341 features are likely attributed to the MFFG algorithm tracking a signal from a sidelobe, and not 
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342 the real wind between these altitudes. Examining the change in wind components over short 

343 time intervals in addition to the wind field itself indicate that ground clutter likely contributed to 

344 the signal around 0100 UTC, 0700 UTC, and 2000 UTC at very low altitudes. Thus, data were 

345 also removed from these regions. The QC’d data are presented in the right panel. 

346 Wind components from LR balloons and the DRWP could also be compared to determine if 

347 the DRWP measurements were acceptable to use on a given day. Balloon data were downloaded 

348 from the KSC Tropical Rainfall Measurement Mission website 

349 (ftp://trmm.ksc.nasa.gov/midds/sonde), and available data were used to visually examine the 

350 characteristics of the wind components from both sources. Following EV44 DOL procedures, 

351 wind components from the closest DRWP profile to 30 minutes after balloon launch were 

352 examined to minimize errors in the comparison associated with the balloon’s rise rate. If the 

353 DRWP profile did not compare well with a balloon profile which was considered acceptable, 

354 then the DRWP measurements would be removed around the time of the comparison. The left 

355 panel of Figure 7 shows that on 11 February 2000 the 1415 UTC DRWP profile deviated from 

356 the 1345 UTC balloon profile above 7 km, with both wind components around 10-15 m s' 1 above 

357 12 km. Thus, DRWP measurements showing this measurement characteristic on this day were 

358 removed. Conversely, the profiles during 11 January 2001 (Figure 7, right panel) show similar 

359 characteristics from both sources. Thus, the DRWP measurements were considered acceptable 

360 on this day. 

361 6. Results 

362 Once the database was QC’d, investigations were performed to examine the algorithm’s 

363 attributes. Missing data tended to exist throughout entire profiles. The ground clutter check, 

364 which is performed manually, mainly affected scattered gates at lower altitudes. However, larger 
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365 


clusters of data contaminated by ground clutter did exist. The vertical beam QC flagged isolated 

366 or narrow vertical regions. The SW, small median, noise power, and isolated datum QC 

367 algorithms all flagged data at sporadic intervals with no general pattern. However, the SW 

368 checks flagged larger clusters of data compared to the other three similar checks. Conversely, 

369 the FGP check flagged adjacent gates at the same altitudes. Gates flagged by the convection 

370 algorithm were removed if data were flagged across extensive vertical regions. The shear checks 

371 tended to flag the boundaries of spurious data regions, with the inside of the regions being 

372 removed manually. The manual QC process also removed vertical discontinuities, sidelobes near 

373 ground clutter and convection, and any other unacceptable feature. 

374 The number and percentage of gates affected by each QC process were tallied. Each process 

375 was assigned a flag, and the number of times an individual flag occurred in each month was 

376 recorded (Table 3). The entire POR contains 162.1 million gates, with a given month containing 

377 12.2 million to 14.3 million gates. Percentages of affected data herein are noted as %POR (% 

378 lowest month to % highest month). Missing data accounted for 35.4% (30.0% to 41.0%) of all 

379 the possible data. The missing data flag was tallied most often because it tended to exist 

380 throughout an entire profile and days in which no data existed were recorded as containing 

381 missing data at every gate and timestamp. The other QC processes combined removed an 

382 additional 6.5% (3.7% to 10.5%) of the possible data. The manual QC process dominated these 

383 QC processes, removing anywhere from 4.8% (1.8% to 8.6%) of the data. The convection QC 

384 process removed 0.6% (0.1% to 1.1%) of the available data. Note that the automated convection 

385 algorithm flagged 2.6% of the data for the POR, but only 0.6% of the data were removed - 

386 indicating the significance of removing flagged data manually. The other automated QC 

387 


processes removed no more than 1.0% of the available data for a given month. The isolated 
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388 datum check removed data at 4,767 gates over the POR. The meteorological wind shear check 

389 removed data at only 570 gates throughout the POR as data that would have been flagged by this 

390 check were likely removed by the DRWP internal shear check. No observations existed that had 

391 unrealistic reports of wind speed or wind direction. The QC’d database contains 58.1% (51.6% 

392 to 64.5%) of the possible wind observations. 

393 Retained complete profiles and pairs were also tallied to support the launch vehicle 

394 community’s interest in examining the vehicle’s entire ascent trajectory. Although the QC’d 

395 database contains profiles which contain data removed by the QC process and can be used for 

396 any application involving winds aloft near KSC, the central focus of generating a DRWP winds 

397 database involved generating a larger sample of complete profiles and profile sets to be used in 

398 vehicle loads and trajectory analyses. 

399 Table 4 depicts the number of complete DRWP profiles and pairs. Generally more profiles 

400 are retained from more recent years than earlier years. March 2000 was the only month over the 

401 POR where zero complete profiles were retained. In addition, the DRWP’s poor performance 

402 during individual periods can be inferred (e.g., February-March, 2000). An average of 30,320 

403 profiles per month exist ranging from 27,436 (October) to 35,239 (March) profiles per month. 

404 No obvious trend in the number of complete profiles seemed to exist from month to month. 

405 Two-hour pairs have an average sample of 15,816 per month ranging from 12,352 (July) to 

406 19,023 (March) pairs per month. Sample sizes for other time separations are on the same order 

407 of magnitude. 

408 The DRWP database has three major advantages over balloon archives. First, the DRWP 

409 database contains on the order of 100 times as many profiles and pairs as the databases derived 

410 from balloon measurements, which would improve confidence in launch simulation results. 
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411 Second, using the DRWP database provides the capability to examine time separations other than 

412 2.0 hours and 3.5 hours as the DRWP pairs archive is not driven by any Program requirements. 

413 In addition, using the DRWP database enables launch vehicle personnel to perfonn assessments 

414 closer than 2.0 hours to launch, which reduces the uncertainty of the wind profile loaded to the 

415 vehicle’s steering commands and the wind through which the vehicle will fly, potentially leading 

416 to launch availability increase due to decreased loads knockdowns. Third, the DRWP database 

417 enables launch vehicle engineers to perform simulations with more than two profiles at a time. 

418 For example; L-3.0 hour, L-1.0 hour, and L-0.0 hour wind triplets can be used to simulate 

419 loading a trajectory at L-3.0 hours, making a GO / NOGO decision at L-1.0 hour, then flying to 

420 the L-0 hour wind. This capability would allow more accurate simulations to be performed 

421 before launch vehicle requirements are written. Times before launch can also be examined to 

422 determine the most optimal DOL assessment sequence. 

423 7. Conclusions 

424 To improve the sample size of MSFC NE’s winds database, QC algorithms were developed 

425 and implemented on DRWP data for the August 1997 to December 2009 POR. A larger sample 

426 of wind measurements not only gives greater confidence in loads and trajectory assessments, but 

427 also provides flexibility to simulate different DOL situations. These features of the DRWP 

428 database should mitigate the limitations of the balloon databases used to support the SSP and 

429 other previous NASA flight vehicle programs. 

430 In addition to increasing the sample size of the database used and providing more flexibility 

431 for DOL simulations in the vehicle design phase, the QC’d DRWP database provides any 

432 upcoming launch vehicle program with the capability to utilize the DRWP profiles on DOL to 

433 compute vehicle steering commands provided the automated and manual QC procedures 
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434 developed are applied to new DRWP data on DOL. In the past, only balloons have been 

435 certified to do this. Although the current DRWP QC process on DOL could be enhanced by an 

436 automated QC process, manual intervention would still be needed to ensure only valid profdes 

437 are used. If high spatial-resolution profiles are desired (such as the SSP’s desire for Jimsphere 

438 measurements) then high frequency components could be randomly added to the DRWP profiles. 

439 The DRWP database provides lots of flexibility in how DOL simulations are performed, and the 

440 QC algorithms provided in this paper will hopefully benefit the aerospace and atmospheric 

441 communities which are interested in utilizing the DRWP. 

442 8. Forward Work 

443 Despite the benefits of utilizing the DRWP database, it does contain a limitation in that only 

444 measurements above 2.7 km are provided for the entire POR. Currently, no QC’d database 

445 exists which contains the sample size of the DRWP database and measurements from near the 

446 surface to 2.7 km. EV44 is thus performing a similar QC to that presented in this paper to data 

447 from the 915-MHz DRWP network at KSC. If an adequate sample exists, complete profiles 

448 from both sources will be combined to generate an extensive database of DRWP profiles 

449 extending from approximately 0.13-18.6 km. The QC algorithms presented in this report, and 

450 any others that are developed for the 915-MHz DRWP, will then be evaluated for operational use 

451 during DOL. 
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Figure 1: Photograph of the KSC 50-MHz DRWP and trailer (left), and a schematic of the 
DRWP’s area and beam configuration (right). 
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Figure 2: 50 MHz Profiler QC 
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Figure 2: GUI used to QC data from the 50-MHz DRWP. A time -height section of the v 
component (m s' 1 ) is plotted in this example. Time is in UTC and altitude is in km. See text for 


descriptions of the different functions. 
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Figure 3: East beam SW (left) and w (right) from 20 August 2001. Time (UTC) is on the x-axis 
and altitude (km) is on the y-axis. SW is colored from 0. 0-2.0 m s' 1 , and w is colored from -2.0- 
2.0 m s' 1 . 
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Figure 4: Convection flag (left) and v (right) on 20 August 2001. Time (UTC) is on the x-axis 
and altitude (km) is on the y-axis. Convection flags are 1 for “convective” and 0 for “non- 
convective”, and v is in m s' 1 . 
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Figure 5: RMS vector error after propagating the FG velocity for each east beam and north beam 
combination. The number of FGPs from the east and north beams are on the x- and y-axis, 
respectively. Colors represent the RMS vector error, and the dashed line is the quadratic fit to 
the black dots, which are maximum points in each FGP combination where the RMS vector is 
less than 1.4142 m s' 1 . 
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Figure 6: Time-height sections of v before (left) and after (right) the QC process is performed. 
Time is on the x-axis in UTC, altitude is on the y-axis in km, and v is in m s’ 1 . 
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Figure 7: Comparisons of u and v from the DRWP and rawinsonde for 1 1 Feb 2000 (left) and 1 1 
Jan 2001 (right). Wind components are on the x-axis in m s' 1 , and altitude is on the y-axis in km. 
DRWP u and v are shown by open circles and squares respectively. Balloon u and v are shown 
by the solid and dashed lines, respectively. 
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575 Tables 

576 

577 Table 1: List of acronyms used. 


DOL 

day of launch 

DRWP 

50-MHz Doppler Radar Wind Profiler 

ESTS 

Engineering, Science, and Technical Services 

EV44 

Natural Environments Branch 

FG 

first guess 

FGP 

first guess propagation 

GUI 

graphical user interface 

KSC 

Kennedy Space Center 

LR 

low resolution 

MATLAB 

Matrix Laboratory 

MFFG 

Median Filter First Guess 

MIDDS 

Meteorological Interactive Data Display System 

MSFC 

Marshall Space Flight Center 

NASA 

National Aeronautics and Space Adminstration 

POR 

period of record 

QC 

quality control 

RMS 

root mean square error 

SNR 

signal-to-noise ratio 

SSP 

Space Shuttle Program 

sw 

spectral width 

u 

zonal wind component 

V 

meridional wind component 
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581 Table 2: Automated QC thresholds. Data were removed if it met the criteria in the threshold 

582 column. 


QC Check 

Threshold 

Vertical SW* 

> 3.0 m/s 

Vertical 0 Doppler Shift* 

| w | > 1.5 m/s and Vertical SNR < 40 dB 

Vertical signal or noise* 

Missing 

Unrealistic wind 

Wind speed <0m/s or wind direction <0 C or wind direction > 360 c 

East or North SW 

> 3.0 m/s 

DRWP Shear 

>0.1 s' 1 

|w | 

> 2.0 m/s 

FGP 

See text 

Me teorlogical Shear 

>0.1 s' 1 

Small Median 

See text 

East or North signal 

Missing 

Rain / Convection 

See text 

isolated datum 

See text 

* Denotes that only data from the vertical beam is re moved 
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Table 3: Number (top) and percentage (bottom) of range gates which were affected by the QC process. Data for each month exists on 
each row and data for each QC process exists on each column. Data matching the criteria in the first three columns were not removed, 
’crcentagcs are rounded to the nearest tenth of a percent. 
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590 Table 4: (Top) Number of complete DRWP profiles assuming at least five minutes between each 

591 profile for each month and year. (Bottom) Number of pairs for each month and time interval, 

592 which ranges from 0.5 hours to 6.0 hours. At least five minutes were skipped between the first 

593 profiles in each pair. 
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